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Abstract 

We investigate form factors of local operators in the multi-component Quantum Non- 
linear Schrodinger model, a prototype theory solvable by the so-called nested Bethe Ansatz. 
We determine the analytic properties of the infinite volume form factors using the coordinate 
Bethe Ansatz solution and we establish a connection with the finite volume matrix elements. 
In the two-component models we derive a set of recursion relations for the "magnonic form 
factors", which are the matrix elements on the nested Bethe Ansatz states. In certain simple 
cases (involving states with only one spin-impurity) we obtain explicit solutions for the 
recursion relations. 

1 Introduction 

One of the goals of many-body quantum physics is the calculation of correlation functions 
of local observables. The form factor program is an approach to tackle this problem; it consists 
of the following three steps: 

1. Finding the eigenstates of the system and inserting a complete set of states between 
the two (or more) local operators. 

2. Evaluating the matrix elements (form factors) of the local operators. 

3. Summing up the resulting spectral series. 

In generic models these tasks present a fantastic challenge. However, the situation is quite 
different in one-dimensional integrable models, where there are exact methods available to 
compute the exact spectrum and the form factors. 

One class of models where this program has been particularly successful is the Bethe 
Ansatz solvable theories related to the sl{2) symmetric i?-matrix, such as the Lieb-Liniger 
model (also known as the Quantum Nonlinear Schrodinger equation) [1, 2] and the XXX 
and XXZ spin chains [3, 4, 5, 6]. A very powerful approach is the Algebraic Bethe Ansatz, 
developed by the Leningrad-school [7, 8, 9], which led to important results concerning the 
scalar products of Bethe states [10, 11] and the form factors [12, 13[. Tremendous effort has 
been devoted to the calculation of correlation functions as well; we do not attempt here to 
review the literature, instead we refer the reader to the book [14[ and the recent paper [15[. 

A different class of models are those non-relativistic theories, where the excitations over 
a fixed reference state have internal degrees of freedom, for example the multi-component 
Non-Linear Schrodinger equation, the sl{N) symmetric spin chains, or the one-dimensional 
Hubbard model. The spectrum of these models can be obtained by the so-called nested Bethe 
Ansatz [16, 17, 18, 19, 20, 21, 22, 23], the algebraic formulation of which was worked out 
for the s/(iV)-related models in the papers [24, 25[. Although the nested Bethe Ansatz is 
successful in finding the spectrum, the construction of the eigenstates is rather complicated 
and there are far fewer results available than in the si {2) case. Norms of eigenstates were 
obtained in [26, 27, 28, 29, 30[ and there are approaches to calculate the scalar products, see 
[31[ and references therein. However, no compact and convenient formulas have been found 
yet, which would facilitate the computation of correlation functions. 
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In the non-relativistic models mentioned above there are exphcit and exact representa- 
tions known for the eigenstates of the system; this ahows (at least in principle) for construc- 
tive methods to find the form factors and correlation functions. The situation is different in 
the realm of (massive) integrable relativistic QFT's [32]. These theories are typically investi- 
gated in infinite volume, the Hilbert-space is spanned by asymptotic scattering states defined 
using the Faddeev-Zamolodchikov algebra [33]. This construction does not allow the direct 
determination of form factors of local operators; an indirect method has been developed in- 
stead: the so-called form factor bootstrap program [34, 35, 36, 37, 38, 39, 40]. The idea is 
to establish the analytic properties of the form factors based on the requirement of locality, 
resulting in a closed set of equations also called "form factor axioms". These equations are 
restrictive enough so that supplied with a few additional assumptions (possibly depending 
on the operator in question) they uniquely determine the form factors. In massive field the- 
ories it is sufficient to obtain explicit expressions for the form factors with a small number 
of particles, because they typically saturate the spectral series for the vacuum correlations 
even at small distances [41, 42]. 

One of the most important form factor axioms is the kinematic pole (or annihilation pole) 
property, which relates (N, M) form factors (matrix elements on an A^-particle and an M- 
particle state) to {N— 1, Af — 1) form factors^. It states that whenever particle rapidities from 
the two states approach each other the form factor has a simple pole (kinematic singularity) 
and the amplitude is given by the form factor with the corresponding particles not present 
and a pre- factor depending on the exact S"- matrix of the theory ]36, 37]. Similar singularity 
properties were also found in the models related to the s/(2)-symmetric i?-matrix in the 
framework of the Algebraic Bethe Ansatz. The pole structure of the scalar products of Bethe 
states was first established in [10], this led to the discovery of the celebrated Slavnov-formula 
]11], a determinant formula describing the scalar product of an eigenstate and an arbitrary 
Bethe state. These results refer to the finite volume states and they were derived using a 
quite general algebraic construction. Moreover, they were used to determine form factors of 
local operators, and in the case when both states are eigenstates, the singularity properties 
of the resulting form factors are found to be essentially the same as in the relativistic case. 
This has been noted recently in [43], where it was also shown that a special non-relativistic 
and small-coupling limit of the sinh-Gordon model form factors yields the known matrix 
elements of the Lieb-Liniger model. 

We wish to note that form factors of local and composite non-local operators in the Lieb- 
Liniger model were also considered using the infinite volume Quantum Inverse Scattering 
Method [44, 45]. One result of this approach is the so-called quantum Rosales expansion, 
which expresses the local field operators using the non-local Faddeev-Zamolodchikov opera- 
tors. The Rosales expansion can be used to read off explicit expressions for the form factors, 
and to establish their analytic properties ]46] , leading to the same kinematical pole equation 
(apart from the normalization) as the one found in [47, 14, 43]. 

The kinematical pole axiom also appears in the study of the form factors of the anti- 
ferromagnetic spin chain ]48]. In this case the states involved are excitations above the 
infinite volume anti-ferromagnetic ground state, which is already filled with a finite density 
of elementary particles. In this respect the situation considered in ]48] is different than in 
]10, 11, 43], where the states in question are elementary excitations over the reference state. 

In this paper we contribute to the calculation of form factors in the multi-component 
Nonlinear Schrodinger equation. Inspired by the results of ]43] we revisit the methods of 
relativistic QFT: we consider the analytic properties of form factors, set up recursion relations 
and make an attempt to find solutions to them, without trying to find manageable expressions 
for scalar products or related, more basic quantities. 

In obtaining explicit solutions to the recursion relations we restrict ourselves to matrix 
elements on states with only a single spin-impurity. One of our motivations to investigate 

^In relativistic field theory the (A'^, M) form factors can be expressed in terms of the analytic continuation of 
the (0, A'^ -I- M) form factors using the so-called crossing relation. Then the kinematic pole is usually expressed in 
terms of the (0, A'^) form factors relating them to the (0, TV — 2) matrix elements. 
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this subclass of form factors is provided by a recent experiment with uhracold atomic gas 
[49], where the motion of spin-impurities was studied in an otherwise polarized background. 
Our present results can form the basis for the theoretical investigation of such situations. 
Related questions were studied in a very recent article [50]. The paper [50] only considered 
the infinite coupling case, whereas our results for the form factors hold at arbitrary coupling 
strengths. 

One of the main steps of the present work is the identification of the (un- normalized) form 
factors in finite and infinite volume. This result bears relevance also for integrable relativistic 
QFT, where related questions have been investigated recently [51, 52]. We give a few remarks 
on this issue in the Conclusions. 

The structure of the paper is as follows. In section 2 we consider the one-component case 
(the Lieb-Liniger model) and establish a number of results about the form factors, using 
only the coordinate Bethe Ansatz wave functions. Although this section does not contain 
new results, it serves as a basis for the generalizations in later sections. In section 3 we recall 
the construction of the (infinite volume) Bethe Ansatz states in the multi-component case and 
we establish the properties of the form factors, in particular the kinematical pole property. 
Section 4 deals with the two-component case: the magnonic form factors are introduced, 
which are the matrix elements on the nested Bethe Ansatz states. A set of "magnonic form 
factor equations" is obtained. In section 5 we solve these equations in a number of simple 
cases, involving states with a single spin-impurity. The sections 3-5 are concerned with the 
infinite volume situation, the connection to the finite volume nested Bethe Ansatz states is 
made in section 6. Finally, section 7 is devoted to our conclusions. 



2 The Lieb-Liniger model: Coordinate Bethe Ansatz 

In this section we review the basic facts about the coordinate Bethe Ansatz solution of 
the one-component Bose gas, the Lieb-Liniger model. The ideas and results of this section 
will be the basis for our investigations of the multi-component systems in section 3, 4 and 6. 



2.1 The model and the coordinate Bethe Ansatz solution 

The second quantized form of the Hamiltonian is 

H= dx (9:^*^9^ ^' + c^'t*!'*^') . (2.1) 

J~L/2 

Here 'i>{x,t) and ^^(^a;,^) are canonical non-relativistic Bose fields satisfying 

[^ix,t),^^y,t)]=S{x-y). (2.2) 
We used the conventions m = 1/2 and h = 1. The first quantized form of the Hamiltonian is 

^ = -Ea^ + 2c^5(a;,-xO- 

The parameter c is the coupling constant; in the present work we only consider the repulsive 
case c > 0. 

In (2.1) L denotes the size of the system. We will consider both the infinite volume 
(L = oo) and finite volume cases. In the latter case we always assume periodic boundary 
conditions. 

The eigenstates of the Hamiltonian (2.1) can be constructed using the Bethe Ansatz 
[1, 2, 14]. The TV-particle coordinate space wave function is given by 



Xn{{x}n 



\{p)n) = 7^ E (-l)'''lexp J z^a;,(7'p), i \{ [(Vp),~{Vp)k-^ce{x,-x^)) , 

■ V&Sn y j J 3>k 

(2.3) 
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where {p} n is the set of pseudo-momenta or rapidities, e(a;) is the sign function, and the 
V G Sn are permutations. The total energy and momentum of the muhi-particle state is 



3 3 

In the infinite volume case the wave function (2.3) is an eigenstate for arbitrary set 
of rapidities. Moreover, the Bethe states with real rapidities form a complete set of states 
[53, 54, 55]. Rapidities with non-zero imaginary parts are not allowed because they result in 
unbounded wave functions. 

In the finite volume case periodic boundary conditions force the quasi-momenta to be 
solutions of the Bethe Ansatz equations: 



^^P,LT^PJ-Pk 



\[ P'-P^-- = l^ j = l...N. (2.4) 



It is known that in the repulsive case (c > 0) considered here all solutions to the Bethe 
equations are given by real rapidities and they provide a complete set of states [56, 57]. 

In finite volume the wave function (2.3) is normalizable. The norm of the eigenstates is 
given by [58, 10] 

A^"(M^)=/ \XN\'^l[{{Pj-Pkr+c')xdetg^^ (2.5) 



j<k 



with 



and 



Gfk = S3,k [L + Y1 'fiiPj -Pl)] - fiPj - Pk) (2.6) 




2.2 Form Factors in finite and infinite volume 

We are interested in the form factors of the field and density operators: 
^-(0) and p(0) = *t(o)$(o). 
The finite volume form factors are defined as 

IP*(MiV,{fc}Ar+i) = 

^ f^/^ (2.7) 

+ 1 / dxi...dxN X*Ni^l:---:^N\{p}N)XN+l{0,Xi,...,XN\{k}N+l) 

J -L/2 



Ki{P}NAk}N) = 

(2.8) 

N dxi. ..dxN-i x^(0,xi, . . . ,a;7v-i|{p}jv)xjv(0, xi, . . . , X7v-i|{fc}jv)- 

J-L/2 

As coordinate space integrals these form factors are well defined for arbitrary rapidities 
and they depend on the volume through the parameters 



— p-VjJ 



The dependence on these parameters was studied thoroughly using Algebraic Bethe Ansatz 
[59, 60, 47, 14[. In particular, it was shown that if the rapidities are solutions to the Bethe 
equations, and there are no coinciding rapidities, then the form factors do not depend on the 
volume explicitly (apart from possible overall phase factors) . 



4 



Note that the form factors are defined using the un-normahzed wave functions, therefore 
the actual finite volume matrix elements can be obtained as 

{{pumku,.) - (2.9) 

and similarly for the density operator. 

An alternative definition for the form factors can be given in infinite volume. In this case 
the real space integrals are not convergent due to the oscillating wave functions. However, 
they can be made convergent by introducing regulators /^(x) in x-space. We choose 

Mx)^fie\x\), 

where f{x) : R"*" — ^ is a smooth function satisfying 

lim f{x) — 1 lim f{x) — 0. 

a;— a:— ^oo 

An example is given by f{x) — e^^ . It can be shown that for every p G K. \ {0} 

lim / dx /e(a;)e*f^ = - lim / dx fe(x)e'P'= = — (2.10) 

independently of the choice of f{x). Actually (2.10) can be considered as a well-defined 
prescription to evaluate the infinite volume integrals. 

Using this prescription the infinite volume form factors are defined as 



/oo ^ ■ 

dxi . ..dxN TT fe{x 



/oo _ _ 

dxi . ..dxN-i Y]_ 



N 

fJx,) 

(2.11) 

X Xjv(2;i, ■ • ■ ,XN\p)XN+liO,Xi, . . .,XN\k) 
N-1 



(2.12) 

X X*Ni^^Xi,. . . ,xn-i\p)xn{0,xi, . . . ,a;Ar„i|fc). 

The connection between the finite volume and infinite volume form factors is established 
by the following theorem. 

Theorem 1. The form factors are the same in finite and infinite volume. In other words, 
if both sets {p} and {k} are solution to the Bethe equations with a volume parameter L and 
there are no coinciding rapidities (pj ^ ki ), then 

Ki{p}N, {k}N+i) = J'^mN, {k}N+l) (2.13) 

Ki{p}N, {k}N) = Ki{p}N, {k}N). (2.14) 

Proof. For simplicity we only consider the field operator and the case = 1 . Then we have 
to prove the equation 

POO 

lim/ dx f,{x) xmP) X2{0:x\ko,k,)^ dx xmp)X2{0,x\ko,ki). (2.15) 



^0 



L/2 



The integrand consists of sums of free wave functions with certain amplitudes. Due to the 
insertion of the field operator the amplitudes depend on the sign of x, therefore the integrals 
have to be split into two parts: 

/-oo /.L/2 /.a /•L/2 

and 



CX3 Jo J-L/2 J-L/2 Jo 
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The integral of an exponential function can be evaluated in the infinite volume case by (2.10) 
whereas in the finite volume case it is given by the Newton-Leibniz formula. The essential 
step to prove (2.15) is to note that in the finite volume case those terms of the Newton-Leibniz 
formula which represent the contributions at x = exactly coincide with the corresponding 
contributions of the infinite volume case. On the other hand, the contributions of the Newton- 
Leibniz formula corresponding to —L/2 and L/2 cancel each other owing to the periodicity 
of the Bethe wave function. Therefore, the two sides of (2.15) are indeed equal. 

Similar arguments can be given in the case of higher particle form factors, and also for 
the form factors of the density operator. □ 

It is clear from this derivation that the case of coinciding rapidities, in particular the 
problem of expectation values has to be treated separately. Whenever there are coinciding 
rapidities the infinite volume FF becomes divergent, the finite volume FF remains finite and 
can be expressed using the properly defined limits of the infinite volume ones [61]. In the 
framework of Algebraic Bethe Ansatz such relations were established for certain non-local 
operators related to correlation functions in [59, 60, 47], whereas local operators describing 
higher-body local correlations were considered recently in [62]. In this work we will only 
consider the case of non-coinciding rapidities. 

2.3 Important properties of the form factors 

The coordinate Bethe Ansatz wave functions (2.3) are completely anti-symmetric with 
respect to an exchange of two rapidities. Therefore 

• ■ • iPAfl^O, • ■ • , + • ■ • , kN) = -J"*(pi, . . . ,PN\kQ, . . . , fcj + i, fcj, . . . , kN) 

J"*(pi, . . . ,pj,pj+i,. . . ,PAr|fco, . . . ,fcAr) = -J"*(pi, . . . ,Pj+i,Pj, . . . ,PAr|fco, . . . , kj^), 

and similarly for the density operator. 

The form factors have kinematical pole singularities whenever pj ki for some j, I. The 
residue of the pole is given by Theorem 2 below. However, before establishing the theorem 
we need the following lemma: 

Lemma 1. Let Dj C M^^^, j ~ 0, . . . , N — 1 be the region 

Xl < X2 < ■ ■ ■ < Xj < < Xj^l < ■ ■ ■ < XJSI-I. 

Then integrating over this region we obtain 

. N-i j . N-l 

Imr / e^i:"--^ H = H U (2-16) 

•^^1 1=1 1=1 l^m=l l=j + l l^m=l f-m 

Proof. The lemma is proven easily by induction, i.e. by performing first the integral over xi 
or a;Ar_i. □ 

Theorem 2. Let p^ fcjv- Then the behaviour of Tfj is given by 
J"*(pi, . . . ,pjv|fco, • ■ • , few) ~ 

^N-l N N-l N \ 

{PN] + ic) WikNj -ic)- ]J {pMj - ic) W{kN] + «c) x (2.17) 

,3 = 1 3 = 1 3 = 1 i=l ) 



kN -PN 



Proof. The form factor is given by a sum over two sets of permutations and a sum over 
the different regions. The pole in p^ — only appears for those permutations when both 
rapidities are coupled to the same Xj . Moreover it follows from Lemma 1 that the singularity 
only appears in those regions, where either Xj is larger than all other coordinates (including 
a; = for the position of the field operator), or if Xj is smaller than all other coordinates. In 
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these cases the integral over xj yields (assuming that the second largest or second smallest 
coordinate is xi) 

lim / dxj /e(a;j)e*('="-P")^^- ^ — ^i{kN-PN)xi 



e^Oj^ " i{kN-pN) i{kN-pN) 



lim / dxj /e(a;j)e'(''"-P")"=^ -> ^ikN-PN)xl 1 



^^OJ-oo iikN-PN) i{kN-pN)' 

Collecting all these terms and adding the proper pre-factors which arise due to the rapidities 
Pn and k^ according to (2.3) we obtain (2.17). □ 

There is an analogous relation for the density operator. The form factors of both operators 
have the structure 

^NmAk}) - n(^^ - ^')n(?" -pj)n^^^ ^ PNmm), (2.18) 

J<1 3>l 3,1 ' 

where ^W({p}|{fc}) are polynomials symmetric in both sets of variables. The degrees of the 
polynomials in their variables are established by the following lemmas. 

Lemma 2. The asymptotic behaviour of the field operator form factor is Pi~^ at pi ^ oo 
and kQ~^ at ko — > oo. 

Proof. At pi — > oo the overall degree of the wave function is p^^^. The leading term factor- 
izes: the amplitude does not depend on the position of the particle with rapidity pi. Therefore 
the integral over the coordinate attached to pi can be performed over the whole real line and 
the regularization scheme yields 

lim / /e(a;i)e'fi^i - lim ( [ f,{x^)e'P'^' + / fe{xi)e'^' A = - + = q. 

Therefore the overall degree of the form factor is determined by the sub-leading terms of 
order p^"^ . The highest possible degree of the coordinate space integrals is pj^^, resulting in 
an overall degree oip^~^. 

At fco ^> oo the leading term in the real space integral for the form factor is given by 
those terms where ko is attached to the coordinate xq — Q. To leading order the reamining 
wave function is proportional to a Bethe Ansatz state with the rapidities {fci, . . . , /cAr} and the 
integrals yield the scalar product (pi, . . . ,pAr|fci, . . . , fc^r). It is assumed that the rapidities are 
different therefore this scalar product vanishes. The sub-leading terms in the wave function 
then yield an overall degree of k^^^ . □ 

It follows from the Lemma above and the factorization (2.18) that in the case of the field 
operator P/y is of degree TV — 1 in all of the pj and the kj variables. Therefore the recursion 
relations (2.17) contain enough information which completely determine the form factors; 
the polynomial Pm can be reconstructed using the Lagrange interpolation procedure [47]. 

In the case of the density operator the polynomial PAr({p}|{fc}), which is symmetric with 
respect to the exchange {p} o {fc}, has the following asymptotic behaviour: 

Lemma 3. The density form factor satisfies the asymptotics 

V FPj{pi...pN\ki...kN) , ..N-i-p^ ( 11 I, \ /oin\ 

Imi^^^ = (-1) Fjv_i(p2 • ■ -PTvlfci • • ■ fcjv)- (2.19) 



Pi 



Proof. The leading terms in the pi — >■ oo limit are given by those permutations where pi 
is attached to the coordinate 0, in other words it is not integrated over. Concerning these 
terms the following relation can be read off from (2.3) 

xn{0,xi,...,xn-i\pi,-.-,Pn) (-1)^"^ , I , , , 

hm j=^XN-i{xi,...,XN-i\p2,---,PN), (2.20) 
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where the sign '-^ indicates that on the Lh.s. only those permutations V £ Sn are kept which 
leave pi at the first place. The statement of the theorem then follows directly from (2.20) 
and the definitions (2.11)-(2.12). □ 

It follows that Pn is of order N in all of its variables. The kinematic recursion relations 
together with the condition (2.19) completely fix the form factor. 

2.4 Determinant formulas for the form factors 

For the sake of completeness we present here explicit determinant formulas for the form 
factors, which will be the basis for the generalizations in section 5. 

In the case of the field operator the form factor can be expressed as [13, 63] 

J'n = Hiikj - hf + c^) dct{S,i - Sn+u). 
Here S* is an x iV matrix defined as 

Y\m={)ik'n ^ k] + *c) Y{m={)0^] ~ ^™ + *'^) 

with 

tiu) — — — 

u{u + ic) 

The first three examples are given explicitly by 

I, , , ^ -^c^{ko-ki){ki-k2){ko-k2){pi-p2) 
■F2 {.Pi,P2\ko,ki,k2j = 



[ko ~Pi){ko -P2){ki -pi){ki ~p2){k2 ~ Pi){k2 - P2) 

X (c^ - koki - kak2 ~ kik2 + (fco + fci + fc2)(pi +P2) - ipiP2)- 

Concerning the density operator, the first determinant formula was established in [12]. 
In the present work we will use an independent representation, which is easily derived from 
the resuhs of [64]: ^ 



TV N 

-I 



m{p}^{k}) = — (-l)~(^+i)/2 [] Y{{ko-pi+ic) X detV. (2.21) 



0=11=1 



Here V is an {N + 1) x {N + 1) matrix with entries 



N 



V, ^ t{k,^p,) + tip,, k,) n _ j,i^i...N 



Vm+i., = T\ 1" ^^t'' V„^+i = l, 

^\ko- Pj + ic 

o— 1 

Vn+i.n+i = 0, 
and 



(2.22) 



i{u) 



u{u + ic) 

We wish to note that (2.21) can be written alternatively as an x determinant, but it is 
useful to keep this form, which makes it possible to find generalizations in section 5. 



^The final form of the determinant formula (2.21) (and the generalizations (5.14) and (5.19)) was suggested 
by Jean-Sebastien Caux. The result given in [64] was expressed as a sum of A'^ determinants, whereas the present 
formula is given by a single determinant, making it more convenient for numerical calculations. 
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The first two cases are given explicitly as 
^?{p\k) = 1 

-^2 Pl'-P2 'Cl,fc2 = 77 771 777 777 7 • 

(fcl -Pl)(fci -P2)[k2 -Pl)(k2 -P2) 

2.5 An alternative representation for the form factors 

The determinant formulas of the previous subsection are very convenient for both an- 
alytical and numerical analysis of the correlation functions. However, it is possible to find 
alternative representations, which might not seem as useful at first sight, but which might 
give clues for the calculation of form factors in nested Bethe Ansatz systems. 

One such representation can be derived using results of the form factor bootstrap pro- 
gram of relativistic integrable QFT's [38]. The papers [65, 66[ considered the form factors 
of breathers (in particular the lowest-lying breathers) in the sine-Gordon model, and they 
arrived at formulas, which give the form factors with a total number of N particles as a sum 
of 2''^ terms. From this result it is possible to derive formulas for the Lieb-Liniger model, first 
performing an analytic continuation in the coupling constant to get the form factors of the 
sinh-Gordon model, and then using a non-relativistic and small-coupling limit as explained 
in [67, 68, 43[. 

In the case of the field operator we obtain the formula 

P T—r 1 

•^*(Miv|{fcWi) = ^ n ^"T^- (2.23) 



Here the polynomial P is given by 



aj=0.iPi=a,l il<i2 

X Y]_ (kii - + (fti - /3j2)*c) W {pil - - (ail - ^2)*^)) 

il<i2 il,i2 



(2.24) 



X 



N N+l 

E(-ir + E(-i)''^ 



The summation is performed over the variables aj = 0, 1 with j = 1 . . . N and l3j = 0, 1, 
with J = . . . iV. In general the number of the aj and l3j variables coincides with the number 
of particles in the bra and ket states, respectively. 

In the case of the density operator the corresponding formula reads 

Ki{p}N\{k}N,f^) = -^^llj^-— (2.25) 

j,k ^ 



with 



P= E E (-1)^^"^+^'* Y[ip^l~P^.2 + {an-a,,)^c)x 

Y[ (kii - + (fti - f3i2)ic) Y[ {Pti - h2 - (a»i - Pi2)icY^ 



X 

il<i2 




(2.26) 



Note that the only difference between formulas (2.24) and (2.26) is the last factor, which is 
called the "p-function" in the original papers [65, 66[ . 
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It can be shown that formulas (2.23)-(2.25) satisfy all necessary conditions established in 
subsection 2.3, therefore they describe the field operator and density form factors indeed. In 
the case of (2.25) the factorization condition (2.19) is also easily checked using 

TV JV \ ^ / ^ ^ 

1 + ^(-1)". + E(-i)''^ + E(-i)"' + E(-i)^^' 

(N N 
^(_l)".+^(_l)ft 
J=2 j=l 



3 Multi-component systems: Coordinate Bethe Ansatz 

In this section we consider the general X-component models in infinite volume. In second 
quantized form the Hamiltonian is 

H= f dx (d.,^]d^-9j + c*|*]*j*/) . (3.1) 

J —oo 

Here 'i>j{x,t) and ^'j(a;,i), j — 1...K are canonical non-relativistic Bose or Fermi fields 
satisfying 

{x, t)^] [y, t) - (y, t)^, {x,t) ^ 5,i5{x ~y), (3.2) 

where a = 1 for bosons (ct = — 1 for fermions), respectively. 
In first quantized form the Hamiltonian is 

^ = -Ea^+2cE'^(^j-^') 

in both cases. Note that the above Hamiltonian is completely spin- independent; the interac- 
tion between the different spin components arises as the effect of the statistics of the wave 
function. 

The construction of the eigenstates of the Hamiltonian (3.1) was established in the papers 
[16, 17, 18, 19, 20] (for a more general scheme see [69]). In the following we collect the main 
results of this procedure; our focus will be on the form factors and their analytic properties. 



3.1 The wave functions 

Before constructing the coordinate Bethe Ansatz wave functions we need to introduce a 
few basic objects and notations. 

Consider the vector space V — . Consider also the iV-fold tensor product 

and a representation p of the permutation group Sn on V^-^^ Let pab denote the action 
corresponding to the elementary exchange {ab). In the physical cases 

Pab — CrPab, 

where Pab is the permutation operator between the vector spaces Va and Vb and a — 1 
(a = —1) in the bosonic (fermionic) case, respectively. 

Consider also a set of parameters {p}n, which will play the role of particle rapidities for 
the Bethe wave function. 

We introduce the operators ]19, 69] 

yab ^ jPj -Pk)Pab - ic 
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Here it is understood that Y°jl^ acts only on the vector spaces a and b and the indices jk 
stand for the rapidities entering the expression (3.3). 

The operators (3.3) satisfy the unitarity condition and the Yang-Baxter equations: 

Y;,' Y,f = 1 (3.4) 

\^ab \^bc -irab _ \^bc \rab -trbc fo r\ 

^jk ^ik ''ij — ^ij '^ik ^jk ■ I'J-Oj 

In the following we attach the rapidities to the vector spaces. To every permutation of 
the rapidities Q e Sn we associate a configuration 

f/'^")^ ® V^'^''^' • • • 4^^^^" . (3.6) 

We define an operator Q(Q, {p}) : (Sn x C^) End(V^(^)) as follows. The permutation 
Q € Sn is re-constructed as a product of elementary permutations and to every exchange of 
rapidities 

we associate the action of Y^'^^; the operator Q{Q, {p}) is defined as the product of the Y^'^'' 
matrices. This definition leads to the property 

Q{Q2, Qip) QiQuP) = Q{Q2Qi,p). (3.8) 

The consistency of the construction is guaranteed by the Yang-Baxter equation (3.5). 
Now we are in a position to construct the vector valued wave functions: 

XN 



The wave functions depend on the (ordered) set of rapidities {p} n and an arbitrary (fixed) 
vector Wat G V^-^^ which is a parameter describing the polarization of the wave function. It 
will be specified in the two-component case in section 4. 
We define the fundamental domain as 

Xj > xi iff j > I. (3.9) 

In this region the wave function is 

Xn{x\p,lun) ^ ^ J2 e'^'^P\-^Q{Q,p)^N. (3.10) 



It can be extended to by symmetry. To write down the formal relation we need the 
representation p{Q) of the permutation Q which is such that 

{Qx)i < ■ ■ ■ < {Qx)n. 
Then the wave function in reads 

XN{x\p.i^N)^^^p{Q-') E e*<''^'«->Q(i?,p)c.^. (3.11) 
^ ReSN 

It can be shown that the wave function defined this way is an eigenstate of the Hamiltonian 
(3.1) for arbitrary {p}n and lon. The total energy and the momentum is given by 



En = ^p] Pn ^J2p3' 



the vector ujn only determines the polarization of the wave function. 
We wish to note that the matrix 

vab _ p vab _ '^(Pj ~ Pk) ~ icPgb 

^]k - ^abljk — (•J-J-^j 

Pj-Pk^' «C 

can be interpreted as the two-particle S-matrix of the theory. Therefore, the individual coef- 
ficients in (3.11) describe the two-particle scattering events. 

It is important to establish the exchange property of the wave function with respect to 
an exchange of rapidities: 
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Lemma 4. 

Xn{x\pi, . . . ,pj,pj+i, . . . ,pn,(^n) = Xn{x\pi, . . . ,pj+i,pj,. ..,PN, Y^j^^ujN). (3.13) 

Proof. It is enough to consider the fundamental region. Let us denote the exchange (j j+1) 
by Pj . We introduce a new summation variable Q' = QPj leading to 



Xjv(a;bi,...,Pj,Pj+i,...,pjv,wjv) = e'^^' ^^p^"'^ Q{Q' Pj,p)wn- 

It follows from (3.8) that 



■ Q'eSN 



QiQ'Pj.p) = Q{Q',PjP)Q{Pj,p) = QiQ',PjP)Y^j^l. 

Therefore 



Xn{x\pi, . . . ,pj+i,pj, . . . ,PN, F^'Z+i^wjv). 



□ 



3.2 The Form Factors 

We are interested in the form factors of the field operators ^P; (0) and the bilinear operators 
Pji{0) = 4'j(0)4';(0). The latter encompass the density operators of particles with a given 
spin (when j = I) and also the spin-flip operators (when j ^ I). 

The (infinite volume) form factors are co-vector valued functions of the rapidities. Eval- 
uated on two vectors ipN € V'^^^ and (j)N+i & they are given as the coordinate space 
integrals 

/■oo 

^Ni{p}N, {k}N+l){lpN, (pN+l) = lim VA^ + 1 / dxi...dxN W fe{Xj)'X 

J-^ f=i (3.14) 



^XAf (xi, . . . , Xn\{p}n, -iPn) C//°VAf+i(0, xi,..., a;jv|{fc}jv+i, 0Ar+i)^ 



N 



N-1 



dxi...dxN-i Y\ f^i^j)^ 

/=i (3.15) 

(uI^\n{p\0,xi,. . . ,xn-i,'>Pn) uPxN{k\0,xi,. . . ^xn-i^^in)) 

\ I N — 1 



Here U^''^ is an operator 



UU) . y(N) ^ yiN-1) (3^^g) 



which acts by taking the scalar product with the unity vector e; in the j-th vector space and 
leaving the others invariant: 

= <5i,Q, (eai (g) ® • • ■ ® eaj_i (g ■ ■ ■ (g ea„) 

The scalar products in (3.14)-(3.15) are the canonical ones in 

and ^^^D, respec- 
tively. In order to conform with our previous notations, in the case of the field operator the 
indexation of the vector spaces corresponding to the rapidities {fcl^v+i is given by 

y(N+i) ^ yko (g,yki (g,...(g, ykN (3^17) 
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In the following we establish the analytic properties of the form factors. The behaviour 
under the exchange of rapidities follows simply from the properties (3.13) of the Bethe Ansatz 
wave functions: 

•^Afbl, ■ • • ,PN\ko, ■ • ■ , kj,kj+i,. . . , kN){(j)N,^N+l) = 

. . . ,P7v|fco, • ■ • , + fcj, • • ■ , fc7v)(07V, S'jj + lV'N+l) 

(3.18) 

J"^(pi, . . . ,p.j,Pj+i,. . . ,PN\ko, . . . , kN){(j)N,tljN+i) = 

F^(pi, . . .,Pj + i,Pj, . . . ,PN\kQ, kN){S^''^^(l)N,1pN+l)- 

Here we introduced the short-hand notation 

Analogous relations hold for the bilinear operators as well. 

The singularity properties of the form factors are more involved. All the poles of the form 
factors arise from the coordinate space integrals and therefore the positions of the poles are 
identical with those in the scalar case, ie. the form factors have poles at pj = ki, whenever 
two rapidities from the two sides coincide. To write down the residues we introduce the 
following notation: 

Idifl® Fn-1 and Fn-i® IdN,N 

are operations where it is understood that a trace is taken with respect to the corresponding 
spaces of V'tv and ipN+i and the remaining vector spaces are substituted in the form factor. 
For example 

(/dl,o ® TN-l){eai ® ■ ■ ■ ® ear,,ebQ® ■ ■ ■ ® e6„) = 

5axM^N-i{(ia^ eQ„,ebj (g) • • • ® eb„) 

With this notation the kinematical poles are given by the following theorem. 
Theorem 3. The pole of the field operator form factor at pn — k^ is given by 

J^'^ipi,. . . ,PN\ko, ■ ■ ■ ,kN){'lpN,4'N+l) ^ -, X 

kN ~Pn 

(J-]v_l (g IdN,N)i4'N, (bN+l) - (T{Idifl ® J-jv_i)(^l'^ . . . 5o,iV . . . Sn-1^N ^N+i) 

(3.20) 

Here we abbreviated = J-'j^_i{pi, . . . ,pN-i\ko, . . . , kpf^i). Poles at other pj = ki can 

be obtained using the exchange property. 

Proof. We follow the ideas of the proof of theorem 2. We consider those terms where pjy 
and k]\[ are coupled to a; at and there is a singularity. Similar to the scalar case, the only two 
possibilities are ii xn is larger, or smaller than any of the other coordinates. To be specific 
we first consider the following two cases: 

xi < ■ ■ ■ < xn^i < xj\j, < Xn and xn < xi < ■ ■ ■ < xn^i, xn < 0. 

In the first case (3.10) gives no action on the two vectors, therefore this term yields 

J-]v-i®Ww,7v. (3.21) 

kN - Pn 

In the second case the wave function is obtained by symmetrization from (3.10): 

Xjv+i(0,xi, . . .,XN\{k}) = Pn-i,n ■ . . PoiXN+iixN,0,xi, . . .,XN-i\{k}). 
Here we used 

(^bi-Pi2 • ■ • Pn~i,n)^ — Pn-i,n ■ ■ ■ Poi- 
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To obtain the coefficient of the exponential 

gi(fcia:iH hfcivaJjv) 

we have to consider the permutation Q e Sm+i 

Q{{ko, kN}) = {fcjv, fco, • ■ • , kN^i} 
For this permutation the corresponding Unear operator is 

Q{Q, {k}) = So.N ■ ■ ■ Sn-i.n- 
Performing similar steps for the dual vector we obtain the scalar product 

a(^PN-l,N ■ ■ ■ -Pl,2<5'^'^ ■ ■ ■ S^~^'^ 4>n\Ui^^ Pn-1.N ■ ■ ■ Pi2Po,1^0.N ■ ■ ■ SN-l.Nft^N+l ^ ^ 

This is equivalent to 

a(^Pl2 ■ ■ ■ Pn-1,nS^''^ ■ ■ ■ ^'^<t)N\Pl2 ■ ■ ■ Pn-1,nUI^^ Sq^n ■ ■ ■ SN-l,N(f>N+l^ ^■ 

Moreover the scalar product is invariant with the permutation of vector spaces within ]/(^) 
therefore the above scalar product is equivalent to 

S^'^ . . . S^^^'^ (IjnIuI'''^ So,N ■ ■ ■ Sn~i,n(Pn+i ^ ^ 

For the contribution in question the above scalar product is equivalent to the action of the 
operator 

(Mi,o <E> T'n-i) (s^ '' ■ ■ ■ S''-^'''cl>N, So,N . . . Sn-i,n<I>n+i) ■ (3.22) 



kN - Pn 

Adding the contributions (3.21)-(3.22), performing the summations over the remaining pos- 
sibilities for the permutations, and using the Yang-Baxter equation one obtains finally the 
statement (3.20). □ 

The residue equation (3.20) involves the parameter a which distinguishes the bosonic and 
fermionic cases. It is useful to write down a relation which does not depend on a. Therefore 
we introduce one more S'-matrix, namely 

Sj,k = (jSj^k S^''' = (jS^'K (3.23) 

Then the kinematical pole equation reads 

^n{P1i ■ ■ ■ iPN\ko, . ■ . , kN){lpN, (pN+i) ^ X 

kN "Pn 

(J-]v_l ® IdN,N)i4'N, (bN+l) - (M,0 ® -F]^_i)(5l'^ . . . So,n . . . Sn-1,N <^Ar+l) 

(3.24) 

It can be shown that an analogous equation (with the S operators involved) holds for the 
bilinear operators and arbitrary higher body local operators as well, irrespective of the statis- 
tics of the model. As a final remark we note that equation (3.24) can be considered as a 
non-relativistic version of the kinematical pole axiom known in the form factor bootstrap in 
integrable relativistic QFT's [36]. 

We conjecture that in the case of the field operator the recursion relation (3.20) together 
with the exchange properties (3.18) determine the form factors uniquely. Then, at least in 
principle they can be constructed with a procedure similar to the one presented in [70] for 
the case of the relativistic 0(3) cr-model. We leave this problem to further research. 

In the case of the bilinear operators it is expected that the recursion relations are not 
restrictive enough. However, similar to the one-component case there is a useful asymptotic 
condition: 
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Theorem 4. The asymptotic behaviour of the form factors of the bilinear operators is given 
by 

lim J^j^{pi, . . . ,pN\ki, . . . ,kN){ipN,4'N) = 

T^_lip2, ■ ■ • ,PN\kl, ■ • • , kN){U^^^1pN, 4>n)- 

Proof. The statement of the theorem is a direct consequence of the relation 

Um '^j{Q)xn{^,Xi,. . . ,XN-l\pi, - ■■ ,PN,i^N) ~ 
Pi— ^oo 



(3.25) 



XN-i(.xi, . . .,XN-i\p2, ■ ■ ■ ,pn,u1^^lon), 



(3.26) 



3 

where similar to (2.20) the sign ^ indicates that on the l.h.s. only those terms are kept where 
pi is attached to xq — 0. Equation (3.26) is checked easily using the definitions (3.10)-(3.11) 
and the the limiting values 



lim Yl'=pab. 



Pj — ^OO 

□ 

We conjecture that the form factors of the bilinear operators are determined uniquely 
by the exchange relations, the kinematical pole axiom (3.24) and the asymptotic condition 
(3.25). Also, we note that (3.25) can be considered as a non-relativistic version of the factor- 
ization property known in relativistic integrable QFT [71]. 

4 Two- component systems: the nested Bethe Ansatz in 
infinite volume 

In this section we consider the two-component models. As a first step we construct the 
so-called nested Bethe Ansatz states. Our approach is somewhat different from the usual one: 
the results presented here apply directly in infinite volume, therefore we are not concerned 
with the periodicity of the wave functions. The connection to the finite volume states is made 
in section 6. 

As a second step we also introduce the "magnonic form factors" which are matrix elements 
of local operators on the nested BA states. We investigate the analytic properties of these 
objects and obtain a set of "magnonic form factor equations". 



4.1 The nested BA states 



We consider the infinite volume Bethe Ansatz states defined in (3.10) in the case of 
V — C'^. The basis in V is formed by the two vectors |-|-) and |— ). Our aim here is to specify 
the vectors ujn entering the Bethe Ansatz states. 

We consider an "auxiliary space" Va = and the operator (also called the "monodromy 
matrix") 

T{u\p) = XaN{u - Pn) ■ ■ ■ Xal(u - pi). (4.1) 

The trace in auxiliary space is called the transfer matrix: 

t{u\p) = TTaT{u\p). (4.2) 

The operator (4.1) can be viewed as the monodromy matrix of an inhomogeneous spin chain, 
where the rapidities pj play the role of inhomogeneities. Then the standard Algebraic Bethe 
Ansatz techniques can be used to construct Bethe states in y . 

To establish the notations we recall that the rational i?-matrix is defined as 



ic 



R{u, c) 



u 
ic 



ic 
u 



\ 



icj 



(4.3) 
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Comparing to the formula (3.12) note that 

mi — ic 

X{u) = Si{u,c)R{u,-cjc) Si{u,c) = — -. (4.4) 

u + ic 

Here Si{u) is the "one-particle" S-matrix, which describes the ampUtude associated to the 
exchange of two particles with the same spin. In the fermionic case it is equal to (— 1), 
whereas in the bosonic case it is just the Lieb-Liniger amplitude. 

In order to conform with the conventions of the spin chain literature we define the nor- 
malized monodromy matrix 

f{u\p) = RaN{u - pn) ■ ■ ■ Rai{u ~ pi). (4.5) 

Here it is understood that in the i?-operators the coupling constant is — crc. With respect to 
the auxiliary space it is written as 



f{u\p) = 



A(u\p) B{u\p) 
C{u\p) D{u\p) 



The commutation relations of the elements of the monodromy matrix can be expressed in 
the compact form 

R{u ~ v)f[u) (g) f{v) = f{v) (g) f{u)R{u - v), (4.6) 

which is understood as an equation in the tensor product of two auxiliary spaces. Eq. (4.6) 
follows from a repeated use of the Yang-Baxter equations (3.5). It follows from (4.6) that 
the transfer matrices t{u\p) form a commuting set of operators. 

We fix the reference state |+)^ — |+ + ■ ■ ■ +)jv ^ T^^^'- The _B(/i)-operators can be 
considered as creation operators of (interacting) spin- waves of |— ) spins acting on the refer- 
ence state. The parameter /i describes the rapidity of the spin wave and is often called the 
magnonic rapidity. 

We define the function ujNi{p}N, {a^Ia/) : x C'*^ ^ V^^^ as 

uJn{{p}n, {m}a/) = B{fii + iac/2\{p}) . . . B{fiM + ^ac/2|{p})|+)^. (4.7) 

The shift of iac/2 is introduced for technical reasons. According to (4.6) the B-operators 
commute with each other, therefore the function a;Ar({p}jv, {/^}a/) is completely symmet- 
ric with respect to the set {/i}. The exchange properties with respect to the set {p} are 
determined once more by the Yang-Baxter relation: 

Lemma 5. The function ll!n({p}n, {j.i}M) satisfies 

un{pi,.-., Pj+1 ,Pj,...,pN,{t^})^ Rj,]+i^N{jPi , • • • , , Vo+1 , ■ • • {a^}), (4.8) 
where Rjj+i = Pj,j+iRj,j+i{pj - pj+i). 

Proof. A modified form of the Yang-Baxter equation (3.5) is 



Y],j+l{P] - Pj+l)Xj+i^a{u - Pj+l)Xja{u - Pj) = 

Xj + i,aiu - Pj)Xjaiu - Pj + l)Yjj+i{pj -Pj+l). 



(4.9) 



This implies 

B{p.\pi, . . . ,pj+i,pj, . . . ,pN)Rj,j+i = Rjj+iB{fi\pi, . . . ,pj,pj+i,. . . ,pn). (4.10) 

By commuting R through the _B-operators and using the fact that R acts trivially on the 
reference state we obtain the statement (4.8). □ 

This leads immediately to the following lemma: 
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Lemma 6. The states defined by (4.7) satisfy the exchange property 

Q{P,p)LUN{{p}N,{f^}M)= n Si{pj-pi)XLUNi{Pp}NAt^}M) (4.11) 



1<1 
Pj>Pl 



for arbitrary P G S*, 



N- 



Proof. It is enough to check the statement for the elementary permutations. Then the state- 
ment follows from equations (4.8) and (4.4). □ 

Finally we obtain the main statement about the nested Bethe Ansatz states: 
Theorem 5. The coordinate space wave functions defined as 

XAr(a;|MjV,MA/)-^ 5] n S,{p,-pi)iON{{Pp}N,{t^}M), (4.12) 

Pj>pi 

where ujn({p}n, {p}m) is given by (4.7), are eigenstates of the Hamiltonian (3.1). Here it 
is understood that (4.12) is defined in the fundamental domain xi < ■■■ < Xm and it is 
extended to the other domains by symmetry. 

Proof. The statement follows from Lemma 6 and the fact that the states (3.10) are eigen- 
states. □ 

It is easy to see that the wave function (4.12) possesses the following exchange property 
with respect to the set {p}: 

Xn{x\pi, . . . ,pj,pj+i,. . . ,PN, {m}a/) = 

Si{pj - pj+i)xN{x\pi,. . .,Pj+i,Pj, ...,PN, {a^Ia/)- 

On the other hand, xn is completely symmetric with respect to the sets {/i}. 

It is useful to obtain explicit representations for the vectors ujN{{p}]sf, {ii}m)- Here we 
just present the known results [47] using the notations of [23]: 

iV 

^n{{p}n,Mm)= ^(ai,-.-,«Af)^-^..ff-"1+), (4.14) 



ai,...,ajvj = l 



where 



A{au...,a,j)^— Y 11 _ {[A{{R^.)ua,l (4.15) 

ReSM l<k<l<M ^ ^ ^'^ 1=1 

where e(a) is the sign function and the propagator of the spin waves is given by 

-i(Tc u - pb + iac/2 

Au, a) = : — - M : — -. (4.16) 

u — Pa — lac/ 2 ^"Ij- u — pb — lacj I 

With this we have finished the explicit construction of the nested Bethe Ansatz states. Note 
that we did not address the completeness of states; in the present approach the magnonic 
rapidities are arbitrary parameters, they are not assumed to satisfy the Bethe equations. 

It is useful to consider the /i oo limit of the wave function. We obtain the following 
statement: 

Lemma 7. Sending a magnonic rapidity to infinity yields the action of the overall spin 
lowering operator: 

lim Hi xn{x\{p}n, ^J'l, ■ ■ ■ , ^J'M) ^ -io-cS- (s^lMw, A^2, • • • , A^a/)- (4-17) 

//I — >oo 
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Proof. It is known from the Algebraic Bethe Ansatz that 

lim II B(ii\{p}n) = -iacS^. (4.18) 

/i— ^oo 

This is easily seen from the construction of the transfer matrix (4.1) or from the explicit 
expression (4.14). Applying (4.18) to (4.7) and finally to (4.12) we obtain the statement of 
the lemma. □ 



4.2 Magnonic Form Factors 

The magnonic Form Factors are defined as the matrix elements of local operators on the 
nested Bethe states. They are scalar functions of four sets of rapidities. 
In the case of the field operator the form factor is defined as 

^n{{p}n, WjM' , {k}N+i, {fijAi) = lim %/7V + 1 / dxi-.-dx^ \\ fe{xj) 

fJi (4.19) 

X (^XNixi, . . .,xn\{p}, {v}) ul"\N+iiO,xi, . . . , XN\{k}, {m}))^- 
Due to spin conservation 

M' = M + 

In the case of the bilinear operators the form factors are defined as 

l-oo N-l 

•^/r(WA'i{i^}M',{fc}Ar,{Ai}Af) = Ihn TV / dxi...dxN~i /eC'z^j) 

j=i (4.20) 

X (yl^^NiO, Xi,..., XN^l\{p}, M) 0nhN{O, Xi,..., XN~l\{k}, {m}))^_^- 

Spin conservation requires 

The Ui''^ operators entering the formulas above are the projectors introduced in (3.16). 

It is important to establish the analytic structure of the magnonic form factors. 

It follows from the properties of the nested wave functions that the form factors are 
completely symmetric with respect to the magnonic rapidities. With respect to the particle 
rapidities they possess the exchange property (4.13) (and its complex conjugate). 

Concerning the kinematical poles we substitute (4.7) into (3.20) and using Lemma 6 we 
obtain 

•^Ar(MI{fc})(wA'(Pl> ■■■,PN, {v}),UJN+l{ko, ■■.,kN, {a*}) 



kN ~Pn 

® IdN.N) (uJNiPl, ■ ■ ■ ,PN, {l^}), il^Ar+l(fco, . . . , /cat, {fJ.})j ~ 
N-l N~l (4.21) 

- n - f"^) n si{PN~pj)x 

3=0 j=l 

cr(/di,0 J'^j^_i)(^UJNiPN,Pl, ■ ■ . ,PN-1,{'^})^ UJN+l{kN, ^0, • • • , k^-l, {/^})) ■ 



X 



Note that in both terms the Id operator acts on those vector spaces to which the rapidities 
piS! and fc^r are attached. Taking the scalar product with respect to these spaces leads to two 
possibilities: Either the corresponding components are + or — . In both cases the resulting 
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amplitudes are evaluated easily using the explicit representation (4.14). We arrive at the 
following "inhomogeneous" form factor recursion equation: 

^N{{p}l..N,W}l..M'\{k}o..N,{lJ-}l..M) ^ -, X 

M M' N-1 N-1 

t=l t=l j=0 k=l 

X ^N-l{{p)l..N-l,{v}l..M'\{k}o..N-lj{lAl--M) 

M' M n 



EE' 



^ (A<r - few - i<7c/2){vs - Pn + icrc/2) 

N-1 N-1 M M' 

n s^{^Jir - k,) n s.{p, - n m^im* - Mr) n ^^^^ - ^*)- 

j=0 j=l t=l t=l 

t^r tjis 

M M' N-1 N-1 

' Il'S'i(Mt - kN)W{flr ~ fJ,t)Y\_^^iPN - yt)W{vt - Us) W Si{kjN) W Si{pNk) 
t=l t=l j=0 fc=l 



X J^N-li{p}l..N-l, {l^}l..s..M'\{k}o..N-l, {^J'}l..f..M)■ 

Here we used 



(4.22) 



^ u - iac u + iac/2 u + iac 

Si(u) — aSi(u) — bi(u) — — W(u) — , 

^ ' ^ ' u + ic 2^ ^ u-iac/2 ^ ' u ' 

and in the last line it is understood that the magnonic rapidities i^s and fi^ are not substituted 
into the form factor. 

Analogous relations can be written down for the bilinear operators as well. In those cases 
the ranges for the different products change according to the number of rapidities present. 

The interpretation of the residue equation (4.22) is the following. The kinematical poles 
arise when two particle rapidities approach each other; in the multi-component case this leads 
to different contributions corresponding to the different spin orientations. In the cases when 
this component is 1+) the remaining part of the wave function has the same number of |— ) 
spins, and its explicit polarization is described by the same sets of magnonic rapidities. This 
corresponds to the second and third lines of (4.22). On the other hand, when the singular 
piece of the wave function carries a |— ) spin, this is associated with one of the magnonic 
rapidities from both states. Therefore the remaining part of the wave functions yields form 
factors with one less number of |— ) spins, ie. one less magnonic rapidity, just as in the form 
factors on the seventh line of (4.22). 

The above recursive equations can be called "inhomogeneous" in the sense that they can 
not be solved by considering fixed sets of {fi} and {i^} as spectator variables: any recursion 
procedure with a given number of magnonic rapidities involves all the form factors with one 
less magnonic rapidity. This makes the solution of the system (4.22) more involved. 

Note that the recursion relation (4.22) has a different structure than the singularity 
properties of the scalar products in the general sl{S) symmetric model considered in [26]. 
This is due to the fact that we considered directly the form factors and not the scalar 
products, and that we used the explicit form of the coordinate wave functions, as opposed 
to the algebraic construction of [26]. Therefore we do not find singularities associated with 
coinciding magnonic rapidities, but the kinematical poles of the particle rapidities yield also 
the "inhomogeneous" terms. 
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The magnonic form factors have the structure 

J'N{{p}N,MM'\{k}N,{^^}M) -C*^'+^^^Jv(MAr,{z^}M'|{fc}jV,{Ai}M)x 

Here Fmin(fc) is the so-called minimal two-particle form factor responsible for the exchange 
properties; it satisfies the relation 

-Fmin(w) = S'i(-u)Fmin(-w)- 

The solutions are given by 

„ , , \ u for fermions 

^^mi„N = <^ „ ^ ^ (4.24) 

— ^ tor bosons. 
^^ u—ic 

Vn is a polynomial which is symmetric with respect to all four sets of variables. 

Lemma 8. The maximal degree oJVn in its variables depends on the operator in question 
and the statistics of the model and is given as follows: 

• In the fermionic case: 

— Field operators: Vn is of order M' in the p variables and of order M — I in the k 
variables. 

— Bilinear operators: Vn is of order M' + 1 in the p variables and of order M + 1 
in the k variables. 

• In the bosonic case: 

— Field operators: Vn is of order N —l-\-M' in the p variables and of order N —l-\-M 
in the k variables. 

— Bilinear operators: Vn is of order N + M' in the p variables and of order N + M 
in the k variables. 

Proof. The total degree of the form factor in the particle rapidities can be established using 
the arguments given in Lemmas 2 and 3 and the explicit form of the wave function (4.12). 
Then the degree of Vn can be read off from (4.23). The main difference between the degrees 
in the fermionic and bosonic cases is a result of the different structure of the minimal two- 
particle form factor. □ 

It follows from the Lemma that in the fermionic case (4.22) is restrictive enough to fix 
the form factors completely. 

On the other hand, in the bosonic case the degree of Vn is typically higher than the 
number of conditions provided by (4.22). Further constraints can be found by sending one 
of the magnonic rapidities to infinity: Lemma (7) and the Wigner-Eckhart theorem provide 
additional relations between different form factors. However, at the present moment it is not 
clear whether (4.22) can be supplemented with other conditions which would fix the form 
factors. These questions are left for further research. 

5 A few explicit solutions for the Form Factors 

In this section we present solutions to the recursive equations (4.22). We only consider 
cases when there is at most one |— ) spin in the two states; for certain operators this leads to 
homogeneous recursive equations which can be solved using generalizations of already known 
techniques. The solution of the full inhomogeneous equations would require new techniques 
and is outside the scope of the present paper. 

As an independent check of our results we also evaluated the coordinate Bethe Ansatz 
expressions for a low number of particles using Mathematica. In the cases iV = 1,2,3 we 



20 



performed the comparisons analytically, whereas in the cases iV = 4, 5 we could only do 
numeric checks. In all cases we found complete agreement. This provides a strong justification 
for the results presented below, especially for the bilinear operators in the bosonic case, where 
the recursion relations do not fix the form factors completely. 



5.1 Fermions 

In the fermionic case we have 5*1 = 1 leading to simple recursive equations. Form factors 
with only |+) spins (no magnonic rapidities) vanish identically, because the one-component 
fermionic model is a free theory (on the level of form factors this follows from S*! = 1 leading 
to vanishing kinematic poles). 



5.1.1 The matrix elements (A^, 0|^' + 1, 1) 

Here we consider matrix elements of the |— ) spin field operator between a completely 
polarized state and a state with only one |— ) spin: 

It follows from (4.22) that the pole at — is given by 

J"jv(pi,...,PAr|fco,...,fcAr,/^) 



kN -PN 

/ _ ^ — kN »c/2 \ jr^ (^p^^ _ _ ^p^_^^i,^^ _ _ (5 1) 



c 



J^N-iiPi^ ■ ■ ■ ,PN-l\ko, ■ ■ ■ , kN-i,n)- 



{kN - PN){kN - M - «c/2) 
The starting point for the recursion is the formal value at = 0: 

^o{ko,ij) = . 

fco — — ic/ Z 

The solution to equation (5.1) is 

-r-( 1 7 ; \ ■ I TT \ ^i<3 ^i>3 



5.1.2 The matrix elements {N, l|^'+|Af + 1, 1) 

If iV > 1 then the "inhomogeneous" term vanishes and we obtain the recursion equations 

J^NiPi'---^PN,'^\ko,---,kN,f^) ^ -, 

kN -PN 

H-kN- ic/2 V-PN + ic/2 \ _^ \i. i, \ /r o^ 

1 1 — ~lo -^n-APi, ■ ■ ■ ,PN-i,'^\ko, ■ ■ ■ ,kN^i, = (5.2) 



fi — kN + ic/2 V — pn — ic/2 



c(pi - v) ^+ 



(fcjv - PAr)(Ai -kN + ic/2){v - Pn - ic/2) 



J^N-iiPi^ ■ ■ ■ ,PN-i,v\ko, ■ ■ ■ , /cjv-i, m). 



At = 1 the inhomogeneous term in (4.22) is the only contribution yielding for the pole at 
ki pi: 

J'i{pi,i^\ko,ki,fi) - 1. , ■ /o —77;{Si{fJ'-ko) - l]-^o^(fco) 

ki — pi ^1 ~ ki + ic/2 u — pi — ic/2 2 

i ic —ic —ic 

ki — pi fi — ki + ic/2 V — pi — ic/2 fi — ko + ic/2 
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There is an analogous relation for the pole at pi — fep- The solution is 
J"i+(pi,i/|fco,fci,/i) = 



ic ic ic fco — ki 



V — pi — ic/2 fi — ki + ic/2 /i — fco + *c/2 (fci — pi)(fco — Pi) 
The solution to the recursive equation then reads 
J^NiPi'--->PN,'^\ko,---,kN,l-J.) = 



c 



fi-kj+ic/2j V - pj ~ ic/2 I Y\, i{kj ~ pi) 



5.1.3 The matrix elements (Af,0|^'^4' |Af, 1) 
The recursive relation is analogous to (5.1): 

J^^^ {pi, ■ . ■ ,PN\ki, . . . ,kN , ij) X 

kN ~Pn 



^N-liPl, ■ ■ ■ ,PN-l\kl, . . . , fcAT-i, /^). 



(5.3) 



{kN - PN){kN - M - ic/2) 
The starting value is 

^t^{pl\kl,^^) = T—r—77;- 

/i — fci + zc/2 

The solution is 

^+-{p„ ...,pN\k„...,kN,f.)^ -^c^ (^{k, - p,) j n ^„fc/+,e/2 fc;^- ^^-'^ 

5.1.4 The matrix elements {N,l\^^_^ ^\N,l) 
The recursion equation reads 

i 



iPi,---,PN,i^\ki,...,kN,tJ-) 



kN -PN 

c{n - v) 



(5.5) 



(fcjv - PAr)(A' - few + ic/2)(i/ - Pn - ■ic/2) 
The starting value is 



-^JV-lbl' ■ • ■ ,PN~l,l^\kl, . ■ . , fcAT-i, ^). 



•^1 (Pi^i^lfciiM) = 



{fi ~ kN + ic/2){i' — pn — ic/2) 
The solution is 

J^N~{Pi, ■ ■ ■ ,PN,'^\ki, ...,kN,fi) = {c{fi - i^))^ ^ 

1 



(5.6) 



det ■ 



HjifJ. - kj + ic/2){v - pj - ic/2) kj - pi 

5.1.5 The matrix elements (iV, l|^^^'+|iV, 1) 

Here the recursion relation is similar to the previous case: 

^n'^{P1, ■ ■ ■ ,PN,l^\kl, . . . ,fcAr,/i) 

r-^jv-ibi' ■ • ■ ^PN-l, v\ki,. . . , kN^l,p) 



(fcjv - Pn){iJ' - kN + ic/2){u -pN - ic/2)' 
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However, this relation is valid only for N > 2. At N = 2 only the inhomogeneous term 
contributes and we have for example the pole at P2 — > ^2 : 

~ T-^ , . ^^^[^1 (M - fcl)-?! (Pl - Z^) - l]^i++(pi|fci) 

k2-p2lJ.-k2 + ic/2y-p2-ic/2^ " ' ^ ^ 

i ic —ic ic{fi — + pi — ki) 

^ P2 fJ- — k2 + ic/2 v — P2 — ic/2 (n — ki + ic/2){v — pi — ic/2) 

There are similar pole relations for the other residues. The solution is 

J^2"^{Pi,P2,'^\ki,k2,n) = 

3 {ki+k2-pi-p2){v- P- + ki + k2-pi-p2) 1 



■det- 



(/i - fci + ic/2){ii - k2+ ic/2){v -pi - ic/2){v - P2 - ic/2) kj — pi 
The solution to the recursive equations is then 



J^N^iPi,- ■ ■ ,Pn, i^\ki, ...,kN,n) = c{c{iJ. - I/)) 



N-2 



X ( - Pj)) ('^ - M + Y^ikj - Pj)) det — ^ 



kj - Pi 



5.2 Bosons 



In the bosonic case Si{v) = (u — ic) / (u + ic) and the recursive equations for the polarized 
states (no magnonic rapidities) produce the form factors of the Lieb-Liniger model. Therefore 
it is expected, that in those cases, where the recursive equations are homogeneous with fixed 
magnonic rapidities, the form factors can be obtained as generalizations of the Lieb-Liniger 
formulas. 

There are three cases where this occurs, with one magnon at most in each of the states: 

• Matrix elements of the down spin field operator: {N, 0\^-\N + 1, 1) 

• Matrix elements of the spin-flip operator: (A'', 0|^'^,_^'_ lAT, 1) 

• Matrix elements of the density of down spins: 

The next matrix elements to consider would be the ones (A/',1|^'^*+|A^, 1). However, 
in this case there is an inhomogeneous term at each step of the recursion; it is given by 
a Lieb-Liniger density form factor. Therefore, already this case is outside the scope of the 
present methods and it requires new ideas. 

5.2.1 The matrix elements (iV, 0|^_|iV -|- 1, 1) 
The residue property at pjv — >■ fcjv+i is 

i 



Tf^{pi,...,PN\ki,...,kN+i,n) , 

kn+1 — Pn 

, N N-l 

/i — k^+i + ic/2 



jx - kN+i - ic/2 ^J-^ 



'Y\_S{kj^N+l) W S{pNk) 



k=l 



-^jV-l(Pl) • • • ,PN-l\kl,.. .,kN,IJ')- 

(5.9) 



The starting point for the recursion is the formal result 

J^Q{ko,n) = IT^^—iK- 

IJL — ko — ic/2 
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It follows from Lemma 8 that the recursion with respect to the p variables completely 
fixes this form factor. We find the following solution: 

where M is an x iV matrix with entries 

Mjk = - MN+i,k 

with 

llm=l "'(.'""1' '''ii llm=l "■('^j' ''^"0 

and 

— r 

h{u) — u + ic hi{u)—u + ic/2 t{u) 



u{u + ic) 
The first example is 

ki — ko ic^^ko + ki — 2/i) 



^1 ip\ko,ki,fJ-) 



ki — ko — ic {fj. — ko — ic/2)(/i — fci — ic/2){p — fco)(p — fci) 



An alternative expression for the same function is a generalization of the formula (2.23): 

•^i^(Miv|{^}w+i,Ai) = 

TT ^ TT ^ TT 1 (5.12) 

cN+1 11 kj -ki-ic J- J- - Pi + ic n, (m - kj - ic/2) i-}kj~pi' 



Here 



P=T. (-1)2:. \{{j>,i-p,, + {a,,^a,,)ic)x 

03=0,1 /3i=0a il<i2 

il<i2 il,i2 

X n(^" + - l)^c/2). 



(5.13) 



It is easy to see that (5.12) satisfies all the conditions and the starting value for the recursion, 
therefore it is equivalent to (5.10). 

5.2.2 The matrix elements (A^, 0|^'^^'_| A^, 1) 

In this case the recursion relation is analogous to (5.9) and we found the following solution, 
which is given by a generalization of (2.21): 



K'{{p}\{k},^i) = l(-i)^(^-i)/2 TT 1 — - TT 1 . 

N yxpsn f.^^j LLk^_ki-tcj^^p,-pi + 2c 

(5-14) 

nn(^o-p,+.c)xdetyx 
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Here V is an (A'' + 1) x (A^ + 1) matrix with entries 
Vji = - M + ic/2)i{kj,pi)+ 



N 



+ {Pi- IJ-- ic/2)i{pi, kj) Yl 



V, 



N 

n 



Po - P] + ic 



0=1 ~ Pi + *^ 
Vn+i,n+i = 



{pi -ko + ic){pi -po- ic) 
{Pi -ko- ic){pi -Po + ic) ■ 

and Vj^N+i = 1, j = 1...N 



j,l = l...N 



(5.15) 



and 



t{u) = 



u{u + ic) 

An alternative representation for the same function is given by 

■^N^i{p}N\{k}N,IJ.) = 



I n 1 

2C^ 11 /; . — — ic 11 p. —p.J^lQ 
J>1 ■' 3>l ■> 



X Px 



-IC 



- kj - ic/2) kj - Pi 



(5.16) 



with 



P= (-1)^^"'+^'^' l[{p^l-p^.+{ai,-a..)ic)x 

aj=0.1 ,a,=0,l il<i2 



il<i2 



il,i2 
N 



(5.17) 



X JJ(^ - fc,- + (2/3,- - l)ic/2) X ^ ((-1)«^- + (-1)'3^) 

The first two cases are given exphcitly as 

T^~{p\k,ii) 



-ic 



jjL — k — ic/2 



J^2 {Pl,P2\kl,k2,n) =(c2 - 4kik2 + 2{ki + k2 -pi -p2)/i+ (fcl +fc2)(pi +P2))x 



k2 - ki 



P2 - Pi 



(/i — fci — ic/2)(/i — k2 — ic/2) k2 — k\ — icp2 — pi+ic 

-c(fci + k2-pi -P2) 



{ki -pi){ki -P2){k2 -Pi){k2 -P2)' 

In the present case the kinematic recursion relation itself is not sufficient to fix the form 
factor. We compared the two representations above to each other and to the coordinate Bethe 
Ansatz results. We performed the analytical check up to = 3 and numerical checks for 
N = 4,5 with Mathematica. Complete agreement was found, which provides a very strong 
justification for the general case of N > 5. 

5.2.3 The matrix elements (A^, 1|*L^'-|A^, 1) 
The residue property at pn — >■ fcjv is 



^jv (pi, • • • ,Piv,z^|fci, . . . , A;jv,/u) 



kN -PN 



1 - 



11 — kN + ic/2 V — Pn — ic/2 
IJ. — kN — ic/2 u — Pn + ic 



N-l 



N-1 



(5.15 



.^jv-i(Pi' • • • >Pjv-i, z^l^i, . . . , kN-i,fJ.)- 
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One solution is 

7-^-(M,.iw,M)=^(-i)-(™nr^^n ' 



c J- ^ kj -ki-ic pj -pi+ic 

J>1 ' 3>l , ^ 

X det F X TT Y[{ko -Pi+ ic) x =— — 

Here V is an {N + 1) x {N + 1) matrix with entries 

Vji = {Pi- + ic/2){jpi - V - ic/2)i{kj,pi)+ 

( ■ i-,\f , • 1 -.-ri {Pi-ko + ic){pi-Po-ic) . 
{pi-li- ic/2){pi -u + ic/2)t{pi,kj) \\ j — ^, 3,1 = 1.. 

VM+i,j = \\ l° ^^t-' Vj,M+i = l, j = l...N 

^Ji ko - Pj + «c 

Visi+1,N+1 = 

and 

t{u) = 



u{u + ic) 
An alternative expression reads 

J^N~i{p}N,l^\{k}N,l^) = 

^-(JV+1) Ej(% -Pj) 



n j,_ _ i-, _ n 



Here 



c2 1 

X P X I I 

n,(M - kj - icl2){v - pj + ic/2) ^.fkj-pi' 



P= E E (-1)^^"^+^'''' l[{Pa-Pi, + {an-ai,)ic)x 

Qj=0,lft=0,l il<i2 

X n ~ ^^'2 + (fti - A2)«c) JJ (Pil - ki^ - (ail - ft2)^c)x 

il<i2 il,i2 



X l[(kj - M - (2/3,- - l)-){pj - V - (2a,- - 1)-) 

i 

The first two cases are given explicitly as 
T^~{p,v\k,\i) = 



c2 



•^2 {Pl,P2,l^\kl,k2,IJ.) 



(5.20) 



1/ - M + 2 (fcj - Pj) f-^i kj -ki-ic J-^J^ Pj -pi+ic (5 21) 



(5.22) 



(5.23) 



{H-k- ic/2){v -p + ic/2) 

c{ki + k2-pi- P2){ki - fc2)(pi - P2) 
{ki -pi){ki -P2){k2 -Pi)(k2 -P2) 

^ 

{11 -ki- ic/2){n -k2- ic/2){v - pi + ic/2)(y -p2 + ic/2) ^ 
{c^{fi - i>) - (fci - k2)'^i> - 2{ki +k2-pi- P2)fii> + +/S(pi - ^2)^)- 
Here we used the auxiliary variables 

ki+k2 - Pi+ P2 

Once again we compared the two representations above to each other and to the coordi- 
nate Bethe Ansatz results up to A'' = 5 and found complete agreement. 
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6 The nested Bethe Ansatz in a finite volume 



The previous sections were concerned with the infinite volume form factors of the multi- 
component systems. Here we consider the two-component case in a finite volume. 

First we recall the conditions for the periodicity of the nested Bethe Ansatz wave function. 

Theorem 6. The wave function (4.12) is periodic in a finite volume L iff the rapidities 
{p}n cind {fJ-}M satisfy the following coupled set of equations: 

^^^^-p^ + u^^^fM_J^^ ; = 1...M. (6.2) 

fJl l^i - P] - «c/2 jj-^ fii - iij - ic 

Proof. For simplicity we only consider the periodicity with respect to the variable xi. By 
symmetry it is enough to consider the periodicity condition 

Xn{-L/2, X2,-.., xn\{p}n, {iAm) ^xn{x2,---, xn, L/2\{p}n, {/^}m) 
in the domain 

< X2 < ■ ■ ■ < Xn < L/2. 
For simplicity we only consider the coefficient of the term proportional to 

gi(2;ipiH hxKPN) g-j 

this will lead to the condition (6.1) with j = 1. At xi — —L/2 the coefficient is equal to 
ajjv({p}, {m})- the other hand, at xi = L/2 the coefficient is 

where Q is the permutation giving Qp = {p2, ■ ■ ■ ,pn,Pi}- Thus we obtain the condition 

^n{{p}, = e'''^'^p{Q-')Q{Q,p)uN{{p}, {fi}). (6.4) 
It follows from the definition (4.1) and from Xa^i{0) = crPa.i that 

T{Pl\{p)) = CrPi2P23 • • • PN-lMPa,NSi,NSi^N-l ■ ■ ■ 81,2- (6.5) 

Here Si,j = y/"^'"' . Taking the trace with respect to the auxiliary space 

t{pi\{p}) = <j''p{Q-^)Si^nSi,n-i ■ ■ ■ Si^2 = <j''p{Q-^)Q{Q,p). (6.6) 
Using (6.4) and (6.6) the periodicity condition is expressed as 

^A^(M, {P}) = e*''^^a^i(pi|M)c.jv(M, {p}). (6.7) 
The eigenvalue equation 

t{u\{p})ujN{{p}, M) = A{u\{p})ujN{{p}, {p}) (6.8) 

can be solved by the standard methods of algebraic Bethe Ansatz. It is known [14] that the 
above equation is satisfied whenever the magnonic rapidities are solutions to the inhomoge- 
neous Bethe equation (6.2). Then the eigenvalues read 

A(.iM) ^ fi ^^(- p.) i n + n n ^ " ^ '^'^'^^ 



ui — u + iac/2 u — p., — iac ui — u + iac/2 

1 ' 1=1 i=\ ' 



(6.9) 
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Substituting u = pi results in 

N M . , 

A(»iM)=n^.(.-p.)n ::::;::: ■ ("»> 

j=i 1=1 ' 

with Si{u) — aSi{u). This leads to (6.1) with j = 1. The other conditions with j = 2, . . . , N 
follow from the symmetry properties of the wave function. It can be checked using the Yang- 
Baxter equation that these equations guarantee the periodicity for the coefficients of all 
exponentials and not only (6.3) considered here. □ 

It is important to note that the states obtained by the finite solutions to the equations 
(6.1)-(6.2) do not span the full Hilbert space: it is known that the Algebraic Bethe Ansatz 
construction (4.7) only produces the states which are highest weight with respect to the 
overall SU{2) symmetry [9]. In order to obtain the non-highest weight states one must act 
with spin lowering operator. This is known to be equivalent (in the proper normalization) to 
sending one magnonic rapidity to infinity. Supplied with these solutions the eqs. (6.1)-(6.2) 
are believed to yield a complete set of states. 

6.1 Form factors in finite volume 

We define the finite volume form factors in the same way as in the case of the Lieb-Liniger 
model. For example, in the case of the field operator the form factor is given by 

rL/2 

^Ni{p}N,{'^}M'Ak}N+lAtJ-}M) ^ VN +1 / dxi...dxN 

J-Lji (6.11) 
X {xn(xx, . . . , xn\{v\, {v^) Ui"\N+iiO, xi,..., XN\{k}, {/i})^ 

Spin conservation requires M' = M + Analo gous definitions can be given for the form 
factors of the bilinear operators. 

Theorem 7. The form factors are the same in finite and infinite volume. In other words, if 
both sets {{plili^}} and {{k}, {fJ-}} are solution to the nested Bethe equations and there are 
no coinciding particle rapidities (pj ^ ki ), then 

I^tvIMjv, {v}M',{k}N+i, {iAm) = T\^{{p}n, {v}M-,{k}N+i, {a*} A/)- (6.12) 
An analogous relation holds for the matrix elements of the bilinear operators. 

Proof. The theorem can be proven with the same arguments as Theorem 1. The key idea is 
that when the wave functions are periodic, the contributions of the Newton-Leibniz formula 
at X = ±L/2 cancel each other and only the contributions at x ~ ±0 remain, which exactly 
coincide with those given by the infinite volume regularization of the oscillating integrals. □ 

Finally we note that the normalized form factors are given by 

{{p}N,MM'\'^l\{k}N+l,MM) - . X f X \ KflSl-X S \ V ^^'^^^ 

V-^(MjV, i^M') TV ({fc}jv+l, {/i}A./) 

and similarly for the bilinear operators. Here A/" denotes the norm of the eigenstates and it 
reads [26, 27, 28] 

M{{p}n. {m}a/) = c'' det g W (l + , , (6.14) 



where Q is an (iV -I- M) x (iV -|- M) matrix, also called the generalized Gaudin-determinant. 
It is the Jacobian associated to the coupled set of nested BA equations and it is given by 



Q = 



_ Vpp Vpi_l 
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where the elements are 



2c' 



1 + cr / ^ v-^ 2c 2c 



2 



7^ C2 + (pj - p;)2 c2 + (pj - pfc)2 

2c' 



(6.15) 



2c \ 2c 



with c' = c/2. 

7 Conclusions and Outlook 

We investigated the form factors of local operators in the multi-component Quantum 
Non-Linear Schrodinger equation. The main results of the present work are the following: 

1. Establishing the analytic properties of the (infinite volume) form factors in the general 
Af-component case; in particular the kinematical pole equation (3.24). This can be 
regarded as a non-relativistic version of the well-known kinematical pole axiom from 
integrable relativistic QFT. 

2. In the two-component case introducing the magnonic form factors and determining 
their analytic structure, in particular the kinematical pole equation (4.22). 

3. The solution of the equation (4.22) in a number of simple cases, involving at most one 
magnonic rapidity per state. This was presented in section 5. 

4. Making a connection to the finite volume form factors, established by the Theorems 1 
and 7, which state that the un-normalized form factors are exactly the same in finite and 
infinite volume. The normalized finite volume matrix elements are given by expressions 
(2.9) and (6.13). 

The most interesting question seems to be whether or not new solutions of the kinematical 
pole equations can be found. A possible direction is to consider integral representations for 
the CO- vector valued form factors in the spirit of the so-called off-shell Bethe Ansatz [37]. 
However, it is not clear if such formulas can be found or if they would be useful for the 
calculation of correlation functions. 

Another interesting question is to consider expectation values of local (or non-local) 
operators in a finite volume. All the techniques presented here apply in the case where there 
are no coinciding particle rapidities; this lies at the heart of identifying the finite volume and 
infinite volume matrix elements. On the other hand, it is known from the one-component 
case that the mean values can always be expressed with the properly regularized diagonal 
limits of the infinite volume form factors. This leads to integral representations for the mean 
values [61, 62]. A natural generalization is to consider this problem in the multi-component 
case; research in this direction is in progress. 

It would be interesting to consider the singularity properties of form factors in the frame- 
work of the Algebraic Bethe Ansatz. To our best knowledge previous works only considered 
scalar products and norms of Bethe states; our kinematical pole equation (4.22) appears to 
be new. A very natural step would be to derive its "finite volume" version from ABA. This 
could help in clarifying the relation between the finite volume expectation values and the 
infinite volume form factors. 

In the end of Section 4 we noted that in the bosonic case the kinematical pole equation 
(4.22) does not contain enough information to determine the form factors. It is an important 
open question whether additional constraints can be found, which would make the recursion 
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equations constraining enough. Also, it is an interesting question why does our formula (5.19) 
work (at least in the cases < 5) for the density operator of the down spin particles. We 
did try other ways to generalize known formulas, and only this one did reproduce the results 
from coordinate Bethe Ansatz. This question is also left for further research. 

As it was already mentioned in the introduction, the infinite volume Quantum Inverse 
Scattering Method (QISM) yields a representation for the field operator in terms of the 
Faddeev-Zamolodchikov operators [44]. From this result (the so-called quantum Rosales ex- 
pansion) the form factors can be simply read off using only the Faddeev-Zamolodchikov 
algebra [46]. The quantum Rosales expansion is established also in the multi-component 
case [72, 73] and this provides us an alternative way to obtain the form factors. We checked 
in a number of simple cases (with low particle number) that the results thus obtained coincide 
with those presented in Section 5 ]74]. However, the QISM does not seem to lead to compact 
and manageable formulas for the generic form factors with higher particle number. On the 
other hand, it would be interesting to consider the two-point functions in this framework: 
it might be possible to derive integral series for the finite temperature and finite density 
correlations along the lines of ]46]. 

In the present work we only considered the non-relativistic multi-component continuum 
systems. However, it is expected that the ideas developed here also apply to other model 
solvable by the nested Bethe Ansatz. In particular it is expected that the annihilation pole 
equation (3.24) (or its "magnonic version" (4.22) and its appropriate generalizations to higher 
rank cases) hold in other sl{N) related models, for example in the SU{N) symmetric Heisen- 
berg spin chains. 

The identification of the finite and infinite volume form factors is an important ingredient 
of the present work. It is very natural to ask: is this result valid in integrable relativistic 
QFT? As it was mentioned in the introduction, in the realm of integrable QFT the infinite 
volume form factors are obtained using the so-called form factor bootstrap program, which 
has been established for theories with both diagonal and non-diagonal S-matrices. On the 
other hand, less is known about the finite volume matrix elements. In [75] it was shown that 
a relation like our (2.9) (concerning the 1-component model) holds in the massive relativistic 
theories with diagonal scattering. However, in the relativistic case the equality is not exact: 
there are finite size effects (decaying exponentially with the volume) due to virtual particle- 
antiparticle pairs "travelling around the world". Concerning non-diagonal scattering theories, 
for example the sine-Gordon model, it is known that the finite size spectrum can be obtained 
(up to exponential corrections) with essentially the same nested Bethe Ansatz construction, 
as applied in this work [51]. The introduction of the "magnonic form factors" is very natural 
also in the relativistic case, and it is expected that a relation like (6.13) holds as well, once 
again up to the exponential corrections. As it was remarked above, an interesting open 
question (both from the coordinate Bethe Ansatz and the integrable QFT point of view) is 
the treatment of the matrix elements with coinciding particle rapidities, or in other words 
the finite volume evaluation of the "disconnected pieces" of the form factors. The resolution 
of this question is important for the evaluation of finite temperature correlations in massive 
integrable QFT ]61, 76, 77, 78]. 
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